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Abstract 

Two,  simple,  solid-oxide  fuel-cell  (SOFC)  power  systems  fed  by  hydrogen  and  methane,  respectively,  are  examined.  While  other  models 
available  in  the  literatures  focus  on  complicated  hybrid  SOFC  and  gas-turbine  (GT)  power  systems,  this  study  focuses  on  simple  SOFC 
power  systems  with  detailed  thermodynamic  modeling  of  the  SOFC.  All  performance-related  parameters  of  the  fuel-cell  such  as  respective 
resistivity  of  the  components,  anode  and  cathode  exchange  current  density,  limiting  current  density,  flow  diffusivity,  etc.  are  all  expressed  as 
a  function  of  temperature,  while  the  flow  through  of  each  nodes  of  the  system  is  described  as  a  function  of  thermodynamic  state.  Full 
analysis  of  the  energy  and  exergy  at  each  node  of  the  system  is  conducted  and  their  respective  values  are  normalized  by  the  lower  heating 
value  (LFIV)  of  the  fuel  and  its  chemical  exergy,  respectively.  Thus,  the  normalized  electrical  energy  outputs  directly  indicate  the  first  law 
and  second  law  efficiencies,  respectively,  of  the  fuel-cell  power  systems.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 

Keywords:  Exergy;  Solid-oxide  fuel-cell;  First  and  second  law  efficiencies;  Fuel  utilization;  Hydrogen;  Methane 


1.  Introduction 

If  one  can  still  recall  the  Kyoto  conference  in  December 
1997  that  highlighted  the  disparate  energy  usage  needs  and 
aspirations  of  many  countries,  he/she  would  have  no  ques¬ 
tion  on  the  imperative  for  clean  and  efficient  use  of  energy 
within  all  countries.  It  is  also  recognized,  in  both  established 
and  fast  growing  economies,  that  the  legitimate  desire  to 
maintain  or  improve  the  quality  of  life  should  be  consistent 
with  the  responsible  use  of  energy.  Responsible  usage  of 
energy  corresponds  to  efficient  use  of  energy  and  minimum 
impact  to  the  global  environment.  To  meet  the  requirements 
of  the  Kyoto  Protocol,  the  fuel-cell  as  an  emerging  technol¬ 
ogy  has  been  considered  to  be  a  potential  candidate  to 
replace  conventional  internal  combustion  engines  for  pro¬ 
pulsion  and  to  integrate  with  gas-turbine  (GT)  technology 
for  power  generation. 

In  power  generation,  because  of  the  synergistic  effects  of 
integrated  solid-oxide  fuel-cell  (SOFC)  and  GT  technolo¬ 
gies,  predicted  results  have  shown  [1,2]  that  an  overall 
system  efficiency  of  70%  (net  ac/lower  heating  value 
(LHV))  or  higher  is  possible  with  a  more  complex  thermo¬ 
dynamic  cycle.  These  studies  were  based  mainly  on  the 
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analysis  of  the  thermodynamic  system  using  first  law  of 
thermodynamics  in  conjunction  with  a  techno-economical 
assessment.  On  other  hand,  increasingly  more  researchers 
are  applying  the  second  law  of  thermodynamics  to  the 
analysis  of  the  overall  plant  efficiency  by  calculating  the 
exergy  content  at  each  node  of  the  thermodynamic  system 
and  the  respective  energy  destruction  in  each  system  com¬ 
ponent  [3,4].  Knowing  the  extent  of  the  losses  distribution  in 
the  thermodynamic  system,  attention  can  be  focused  on 
improvements  in  a  particular  system  component  or  process. 
A  detailed  description  of  exergy  analysis  and  the  reduction 
of  irreversibility  has  been  given  by  Cornelissen  [5].  This 
study  also  provided  a  basis  for  the  relationship  between 
exergy  analysis  and  environmental  life-cycle  analysis 
(LCA). 

In  this  paper,  we  present  two  simple  SOFC  power  systems 
fed  by  hydrogen  (FF)  and  methane  (CH4),  respectively.  A 
simulation  technique  is  applied  and  individual  system  com¬ 
ponent  models  are  developed  and  integrated  in  the  simula¬ 
tion  code.  Emphasis  is  placed  on  the  development  of  a  more 
complete  SOFC  model,  which  can  be  used  to  simulate 
different  designs  of  SOFC  provided  that  the  geometric 
specifications,  performance-related  parameters  and  operat¬ 
ing  conditions  are  specified.  The  study  also  provides  the 
foundation  for  future  analysis  of  a  more  complex  power 
system. 
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Nomenclature 

a  coefficient 

C  concentration 

e  energy  of  molecular  interaction 

E  voltage 

Ex  exergy 

F  Faraday  constant 

g  acceleration  due  to  gravity 

H  enthalpy 

i  current  density 

ia  exchange  current  density 

/oa  anode  exchange  current  density 

/oc  cathode  exchange  current  density 

I  current 

I  distance 

/a  distance  from  anode  surface  to  reaction  sites 

Zc  distance  from  cathode  surface  to  reaction  sites 

M  molecular  mass 

n  mole  number 

p  pressure 

Q  heat 

R  universal  gas  constant 

S  entropy 

T  temperature 

V  velocity 

W  work  done 

Greek  letters 
t]  polarization 

tj Act  activation  polarization 

tjconc  concentration  polarization 

i/ohm  Ohmic  polarization 

p  resistivity 

a  collision  diameter 

ercv  entropy  production 

O,)  collision  integral  based  on  the  Lennard-Jones 

potential 
c  tortuosity 


2.  Exergy 

2.1.  Exergy  concept 

Exergy  is  defined  as  the  maximum  obtainable  work  a 
substance  can  yield  when  it  is  brought  reversibly  to  equili¬ 
brium  with  the  environment.  The  exergy  method  of  analysis 
is  essentially  based  on  the  second  law  of  thermodynamics 
and  the  concept  of  irreversible  entropy  production.  It  is 
concerned  with  the  efficiency  of  use  of  the  available  work 
that  is  generated  from  energy  resources,  and  therefore, 
enables  determination  of  the  location,  type  and  true  magni¬ 
tude  of  each  loss  in  a  thermodynamic  system.  Such  informa¬ 
tion  is  useful  when  designing  a  thermal  system  or  reducing 
sources  of  inefficiency  in  an  existing  system. 


The  environment  is  generally  the  surrounding  environ¬ 
ment  of  the  system  since  the  final  deposition  of  the  initially 
available  energy  will  be  through  external  cooling  devices  or 
ejection  of  an  energy-bearing  substance  into  the  environ¬ 
ment.  It  is  a  very  large  reservoir  in  the  state  of  perfect  ther¬ 
modynamic  equilibrium,  in  which  no  gradients  of  pressure, 
temperature,  chemical  potential,  kinetic  or  potential  energy 
can  exist.  Thus,  there  is  no  possibility  of  producing  work 
from  any  form  of  interaction  between  parts  of  the  environ¬ 
ment.  The  environment,  therefore,  is  a  natural  reference  for 
assessing  the  work  potential  of  different  kinds  of  systems. 

Three  forms  of  exergy  transfer  are  usually  established  to 
perform  an  exergy  analysis  of  a  system,  namely,  exergy 
transfer  with  work  interaction,  heat  interaction  and  material 
streams.  Other  forms  of  exergy  transfer  include  friction, 
momentum,  potential  interaction,  etc. 

In  conducting  exergy  analysis,  two  forms  of  equilibrium, 
the  environmental  state  and  the  dead  state  are  considered. 
The  environmental  state  is  a  restricted  equilibrium  where  the 
conditions  of  mechanical  (p)  and  thermal  ( T)  equilibrium  are 
satisfied.  The  dead  state  is  an  unrestricted  equilibrium  where 
the  conditions  of  mechanical  (p),  thermal  (7)  and  chemical 
(/<)  equilibrium  are  satisfied.  Under  this  full  thermodynamic 
equilibrium  between  the  system  and  the  environment,  the 
system  cannot  undergo  any  changes  of  state  through  any 
form  of  interaction  with  the  environment. 

The  control  volume  of  a  thermodynamic  system  with  a 
single  inlet  and  outlets  is  shown  in  Fig.  1 .  The  exergy  of  this 
system  can  be  express  as: 


E'-P 


+  m 

CTcv  =  0 


Wc 


(hi  —  hi)  —  To  (Si  —  S2)  H — — ^ — ~  +  g(zt  —  Z2) 


(1) 


For  each  component  process,  an  exergy  balance  can  be 
established.  By  considering  all  the  in-going  and  out-going 
exergy  flows,  the  exergy  destruction  can  be  calculated. 
When  dealing  with  the  exergy  of  a  process  component, 
the  difference  between  exergy  losses  and  destruction  should 
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Fig.  1.  An  open  thermodynamic  system  with  a  single  inlet  and  outlet. 
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be  noted.  Exergy  losses  consist  of  exergy  flowing  to  the 
surroundings,  whereas,  exergy  destruction  indicates  the  loss 
of  exergy  within  the  system  boundary  due  to  irreversibility. 

2.2.  Exergy  components 

In  calculating  the  exergy  of  a  stream  of  matter,  the  exergy 
can  be  divided  distinctively  into  different  components.  In  the 
absence  of  nuclear  effects,  magnetism,  electricity  and  sur¬ 
face  tension: 

Ex  =  Ek  +  Ep  +  Ephy  +  £chem  (2) 

where  Ek,  £p,  Ephy,  and  £chem  are  kinetic  exergy,  potential 
exergy,  physical  exergy  and  chemical  exergy,  respectively. 

In  this  study,  only  physical  and  chemical  exergies  are 
considered  at  different  nodes  of  the  system,  while  kinetic 
and  potential  effects  are  neglected.  Physical  exergy  is 
defined  as  the  maximum  amount  of  work  obtainable  when 
the  stream  of  substance  is  brought  from  its  initial  state  to  the 
environmental  state  defined  by  p0  and  T0,  by  physical 
processes  which  involve  only  thermal  interaction  with  the 
environment.  Chemical  exergy  is  defined  as  the  maximum 
amount  of  work  obtainable  when  the  substance  is  brought 
from  the  environmental  state  to  the  dead  state  by  processes 
involving  heat  transfer  and  exchange  of  substance  only  with 
the  environment. 


3.  Modeling 

The  simulation  program  developed  consists  of  two 
mathematical  models,  namely,  the  electrochemical  model 
and  the  heat  exchanger  model  for  fuel-cell  systems  operat¬ 
ing  on  hydrogen  or  methane.  A  PE-fed  SOFC  system, 
which  incorporates  two  pre-heaters,  a  SOFC  stack  and  an 
afterburner,  is  shown  in  Fig.  2.  The  slightly  more  complex 


CPE-fed  SOFC  system,  as  shown  in  Fig.  3,  incorporates  a 
mixer,  a  vaporizer,  two  pre-heaters,  an  external  reformer, 
a  SOFC  stack,  and  an  afterburner.  The  reformer  used  is  a 
steam-reforming  type,  which  requires  an  external  heat 
source  for  the  endothermic  reaction.  The  reformer  reforms 
the  CH4  to  PE-rich  reformate  following  the  chemical 
reaction: 

CH4  +  2.444PEO  -»•  C02  +  4H2  +  0.444H,O  (3) 

The  minimum  mole  number  of  water  vapor  required  for 
steam  reforming  is  two.  The  additional  0.444  moles  of  water 
vapor  are  used  to  prevent  possible  thermal  decomposition  of 
the  methane  that  would  otherwise  form  solid  carbon  which 
would  deactivate  the  catalytic  reaction  of  the  reformer. 

The  simulation  program,  written  in  Visual  Basic,  allows 
the  user  to  specify  the  operating  pressures,  fuel  supplying 
rate,  fuel  utilization  rate,  percentage  of  theoretical  airflow, 
reformer  efficiency  and  temperatures  at  certain  points  of  the 
system.  After  accepting  the  user’s  entry  data,  it  performs  a 
calculation  to  determine  the  molar  chemical  compositions  of 
the  flow  streams  at  different  nodes  of  the  system  and  the 
amount  of  air  required  by  the  SOFC  stack  and  afterburner  to 
sustain  the  system  operation.  Then,  it  calls  on  the  electro¬ 
chemical  model  to  determine  the  electrical  work  output  of 
the  SOFC  stack,  its  associated  heat  rejection  to  the  surround¬ 
ings,  and  the  exit  flow  temperature.  Following  this,  the  heat 
exchanger  model  is  called  to  determine  the  exit  temperatures 
of  the  two  pre-heaters.  Based  on  these  calculated  tempera¬ 
tures,  the  enthalpy  and  exergy  at  corresponding  points  can  be 
determined.  These  enthalpies  and  exergies  are  then  normal¬ 
ized  relative  to  the  FHV  and  chemical  exergy  of  the  fuel, 
respectively. 

The  physical  and  chemical  exergies  are  calculated  using 
the  following  equation: 

Fphy  =  £chem  =  (Hn  -  H°)  -  T0(S„  -  S°)  (4) 


To 


Fig.  2.  Hydrogen-fed  fuel-cell  system. 
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Water 


Fig.  3.  Methane-fed  fuel-cell  system. 


where  subscript  n  refers  to  the  state  points  of  a  system, 
°  refers  to  environmental  condition,  H  is  the  enthalpy,  and 
S  the  entropy. 

Assuming  that  the  gases  obey  ideal  gas  behavior,  their 
respective  enthalpies  and  entropies  can  be  determined  from 
the  following  polynomial  equations  fitted  to  the  data  of  the 
JANAF  table  [6], 


^  =  ai+Tr+fr2+Tr5+yr*+?  (5) 

T  =  0|ln(r)  +a2r  +  ‘jT2  +  ‘jT3  +  ‘jT4  +a7  (6) 


The  coefficients,  cq-t/y,  for  different  gases  used  in  this  study 
can  be  found  in  Table  1. 

The  assumptions  and  conditions  of  model  used  in  the 
simulation  program  are  listed  as  follows: 


•  steady  flow  with  negligible  frictional  losses; 

•  negligible  changes  of  potential  and  kinetic  energies  in  any 
processes; 

•  the  environment  is  at  STP  conditions,  i.e.  298  K  and  1  atm; 

•  isothermal  process  in  the  afterburner; 

•  all  other  gases  except  hydrogen  and  oxygen  are  consid¬ 
ered  to  be  diluents; 

•  wet  atmospheric  air  is  used; 

•  stoichiometric  fuel-air  reaction  throughout; 

•  default  fuel  supplying  rate  is  1  kg  s~ 1 . 


3.1.  Electrochemical  model 


The  electrochemical  model  considers  the  effect  of  all 
forms  of  overpotential  on  the  SOFC  performance,  its  asso¬ 
ciated  heat  lost,  and  the  electrical  energy  produced.  Since 
all  the  performance  related  parameters  are  dependent  on 


temperature,  iterations  are  performed  to  calculate  the  stack 
temperature,  with  known  thermodynamic  state  at  the  inlets 
of  the  fuel-cell,  until  convergence  criteria  is  met. 


3.1.1.  Activation  overpotential 

To  avoid  the  ambiguity  of  simplified  models,  such  as  the 
Tafel  equation  or  a  linear  potential-current  relation,  used 
under  different  operating  conditions,  the  following  general 
Butler- Volmer  equation  is  used  to  calculate  the  respective 
overpotential  of  the  anode  and  cathode: 


i  =  ;'0<J  exp  (  0  '7e^;v'  )  -  exp 


-(1-/0 


»eT»/Act' 


RT 


(7) 


where  ft  is  the  transfer  coefficient  and  iQ  the  ‘apparent’ 
exchange  current  density.  The  transfer  coefficient  is  con¬ 
sidered  to  be  the  fraction  of  the  change  in  polarization  that 
leads  to  a  change  in  the  reaction  rate  constant;  its  value  is 
usually  0.5  in  the  context  of  a  fuel-cell.  Hence, 


i  =  2  i0  sinh 


f  "eFii  Ac t\ 

V  2 RT  ) 


(8) 


or 


2  RT 

'/Act  =  —p  Slnh 

neF 


(9) 


The  ‘apparent’  exchange  current  density  is  a  function  of 
temperature,  which  is  related  to  the  charge  transfer  resis¬ 
tance  by: 


At  = 


RT 

nFi0 


(10) 


or 

RT 

L  = - 

nFr ct 


(ID 


Table  1 

Coefficients  used  in  Eqs.  (5)  and  (6)  for  various  gases3 


Species 

Coefficients 

a\ 

d'l 

CI4 

a6 

df 

300  K  <T<  1000  K 

Ar 

2.50000000E+00 

0.00000000E+00 

0.00000000E+00 

0.00000000E+00 

0.00000000E+00 

— 7.45375020E+02 

4.36600060E+00 

ch4 

1.50270720E+00 

1.04167980E— 02 

—  3.91815220E— 06 

6.77778990E— 10 

— 4.42837060E— 14 

— 9.97870780E+03 

1.07071430E+01 

C02 

4. 460804 10E+00 

3.098 17 190E-03 

— 1.239257 10E— 06 

2.27413250E— 10 

— 1.55259540E— 14 

— 4.89614420E+04 

— 9.86359820E— 01 

h2 

3.10019010E+00 

5.1 1 194640E— 04 

5.26442100E— 08 

— 3.49099730E—  1 1 

3.69453450E— 15 

—  8.77380420E+02 

— 1.962942 10E+00 

H20  (G) 

2.71676330E+00 

2.945 13740E-03 

—  8.02243740E— 07 

1.02266820E— 10 

— 4.84721450E— 15 

— 2.99058260E+04 

6.63056710E+00 

h2o  (L) 

1.27127820E+01 

—  1 .76627900E— 02 

— 2.25566610E— 05 

2.08209080E— 07 

— 2.40786140E— 10 

— 3.74832000E+04 

— 5.91153450E+01 

n2 

2.8963 1940E+00 

1.51548660E— 03 

— 5.72352770E— 07 

9.98073930E— 11 

— 6.52235550E— 15 

— 9.05861840E+02 

6.16151480E+00 

02 

3.62195350E+00 

7.361 82640E-04 

—  1 .96522280E— 07 

3.620155880E— 11 

— 2.89456270E— 15 

—  1.20198250E+03 

3.61509600E+00 

5000  K  <  T  <  1000  K 

Ar 

2.50000000E+00 

0.00000000E+00 

0.00000000E+00 

0.00000000E+00 

0.00000000E+00 

— 7.45374980E+02 

4.36600060E+00 

CH4 

3.82619320E+00 

— 3.97945810E— 03 

2.45583400E— 05 

— 2.27329260E— 08 

6.96269570E— 12 

—  1.01449500E+04 

8.66900730E— 01 

co2 

2.40077970E+00 

8.73509570E— 03 

— 6.60708780E— 06 

2.00218610E— 09 

6.32740390E— 16 

— 4.83775270E+04 

9.695 14570E+00 

h2 

3.057445 10E+00 

2.67652000E— 03 

— 5.80991620E— 06 

5.52103910E— 09 

—  1.81227390E— 12 

— 9.88904740E+02 

— 2.29970560E+00 

H20  (G) 

4.07012750E+00 

—  1.1 0844990E— 03 

4.15211800E— 06 

— 2.96374040E— 09 

8.07021030E— 13 

— 3.02797220E+04 

— 3.22700460E— 01 

h2o  (L) 

1.27127820E+01 

—  1 .76627900E— 02 

— 2.25566610E— 05 

2.08209080E— 07 

— 2.40786140E— 10 

— 3.74832000E+04 

— 5.91 153450E+01 

n2 

3. 674826 10E+00 

— 1.208 15000E— 03 

2.32401020E— 06 

— 6.32175590E— 10 

— 2.25772530E— 13 

- 1.061 15880E+03 

2.35804240E+00 

02 

3.62559850E+00 

— 1.8782 1840E— 03 

7.05545440E— 06 

-6.7635 1370E-09 

2.15559930E— 12 

-1.04752260E+03 

4.30527780E+00 

a  JANAF  table. 
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where 

io=Ae~WRT)  (12) 


3.1.2.  Ohmic  overpotential 

Ohmic  overpotential,  which  contributes  by  the  electro¬ 
lyte,  electrodes  and  interconnector  of  the  fuel-cell,  occurs 
because  of  the  resistance  to  the  flow  of  ions  in  the  ionic 
conductors  and  the  resistance  to  electrons  through  the 
electronic  conductors.  Since  these  resistances  obey  Ohm’s 
law,  the  overall  Ohmic  overpotential  can  be  written  as: 


*7  Ohm  —  iR 


(13) 


The  resistance  of  each  material  used  in  the  SOFC  com¬ 
ponents  can  be  calculated  from  its  respective  resistivity, 
which  is  a  function  of  temperature  [7],  i.e. 


and 


(14) 


p  =  ae ^ 


(15) 


where  a  and  b  are  the  constants  which  are  specific  to  the 
materials. 


3.1.3.  Concentration  overpotential 

To  avoid  a  heavily  reliance  on  correlation  to  determine  the 
limiting  current  density,  the  complete  concentration  over¬ 
potential  equations  derived  in  a  previous  study  can  be  used 
[8],  namely: 

I  -  (RT/2F)(ld/Dcff/dp]H2)i' 
1  +  {RT/2F)(ljDeftaplH20)i 
(16) 


RT 

Jonc,a  —  2^7 


Knudsen  diffusion,  molecules  collide  more  frequently  with 
the  pore  walls  than  with  other  molecules.  Upon  collision,  the 
atoms  are  instantly  adsorbed  on  to  the  surface  and  are  then 
desorbed  in  a  diffused  manner.  As  a  result  of  frequent 
collisions  with  the  wall  of  the  pore,  the  transport  of  the 
molecules  is  impeded. 

In  the  Knudsen  diffusion,  the  diffusion  coefficient  for 
straight  and  round  pores  is  expressed  by  [9]: 


To  account  for  the  tortuous  path  of  the  gas  molecules 
rather  than  travel  along  the  radial  direction,  and  for  the 
porosity  of  the  ceramic  given  that  diffusion  occurs  only  in 
the  pore  space,  an  effective  Knudsen  diffusion  coefficient  is 
defined: 


Ayeff  =  £>KT  (19) 

£ 

Using  the  Chapman-Enskog  theory  of  prediction,  the 
binary  ordinary  diffusion  coefficient  in  the  gas  phase  can 
be  calculated  by  [10]: 


Db  =  1.8583  x  10“'  —  +  — 


1 


1  \  L/2  p3/2 


Ma  Mb/  /wab--d.ab 


(20) 


where  p  is  the  total  pressure  in  atm;  rrAB  =  (erA  +  erB)/2  the 
collision  diameter  in  angstroms  (A);  QD  AB  the  collision 
integral  based  on  the  Lennard-Jones  potential  and  can  be 
obtained  from  eAB  which  is  the  energy  of  molecular  inter¬ 
action  in  ergs  (eAB/k  =  sj  {eA/k)(eB/k)). 

Similar  to  Knudsen  diffusion,  the  effective  diffusion 
coefficient  for  ordinary  diffusion  is  defined  as: 

AxB.eff  =  -Dab  t  (21) 

£ 


cathode  :  rjConc  c 


RT  ^  {pc/do2  -  {pc/So2  -  pp2)  exp(RT/4F(8o2lc/Deff}Cpc)i)) 

2 F  n  P\h 


(17) 


Concentration  overpotential  becomes  significant  when 
large  amounts  of  current  are  drawn  from  the  fuel-cell. 
The  partial  pressures  of  the  gases  at  the  reaction  sites,  which 
corresponds  to  the  volume  concentrations  of  the  gases,  will 
be  less  than  that  in  the  bulk  of  the  gas  stream  when  a  large 
amount  of  current  is  drawn.  If  such  partial  pressures  or 
concentrations  of  gases  are  unsustainable,  concentration 
polarization  will  cause  excessive  voltage  losses  and  the 
fuel-cell  will  cease  to  operate. 

Diffusion  through  a  porous  media  can  be  described  by 
either  ordinary  or  Knudsen  diffusion.  In  this  context,  both 
types  of  diffusion  are  included  in  the  model  to  render  it  more 
generalized.  Ordinary  diffusion  occurs  when  the  pore  dia¬ 
meter  of  the  porous  material  is  large  in  comparison  with  the 
mean  free  path  of  the  gas  molecules,  whereas,  for  Knudsen 
diffusion,  molecular  transport  is  through  pores  which  are 
small  in  comparison  to  the  mean  free  path  of  the  gas.  In 


To  account  for  the  combined  effect  of  Knudsen  and 
ordinary  diffusion  in  the  diffusion  process,  a  combined 
diffusion  coefficient  is  used,  i.e. 


1  _  i  1 

DAjeff  U,\B.eff  DAK.eff 


(22) 


3.1.4.  Electrical  power  and  heat  transfer  losses 

Once  all  the  overpotentials  are  computed,  the  cell  voltage 
can  be  determined  by  subtracting  these  overpotentials  from 
the  Nernst  potential: 

T(f  )  To  (*/ohm4~  */Act,a4~  */Act,c"6  *7Conc,a4~  */Conc,c)  (23) 

where  the  Nernst  potential  is  given  by: 


RT  RT 

=  —  In  K - In 

2  F  4  F 


(24) 
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Table  2 

Conditions  of  heat  exchanger11 


Type  of  fuel 

Pre-heater 

no. 

Inlet  temperature 
cold  side  (K) 

Outlet  temperature 
cold  side  (K) 

Inlet  temperature 
hot  side  (K) 

Outlet  temperature 
hot  side  (K) 

Calculated  pre-heater  size 
(m2)  (include  fouling  factor) 

Hydrogen 

2 

298 

1100 

1200 

1157 

650 

1 

298 

1117 

1157 

382 

63400 

Methane 

2 

298 

1140 

1200 

1143 

550 

1 

298 

1102 

1143 

412 

26300 

11  Operating  conditions:  fuel  flow  rate,  1  kg  s  air  flow  rate,  600%;  theoretical  air,  1  atm  (both  anode  and  cathode);  fuel  utilization  rate,  75%;  reformer 
efficiency,  90%  (methane-fed  system  only);  overall  heat  transfer  coefficient,  0.05  kW  m~2  K. 


The  electrical  power  density  produced  by  the  fuel-cell  is 
calculated  by: 


W  =  El 


(25) 


capacity  rates  is  then  available.  Thus, 
Umin 


Cr  = 


C„ 


(31) 


Similarly,  the  heat  lost  to  the  environment  at  an  equili¬ 
brium  state  can  be  determined  readily  by  evaluating  the 
entropy  rate  balance  for  a  control  volume,  i.e. 

Qcv  =  T(AS  -  ercv)  (26) 

where 


AS 


(27) 


The  ‘thermodynamic’  entropy  production,  which  is  actu¬ 
ally  the  irreversibility,  is  related  to  the  ‘electrochemical’ 
overpotential  as: 


The  effectiveness  of  cross-flow  heat  exchangers  can  be 
calculated  from: 


1  _  e[— NTU(1-Cr)] 

1  -  Qel-NT^1-^)]  ’ 


forCr  <  1 


NTU 

1  +  NTU  ’ 


for  Cr  =  1 


where 


(32) 

(33) 


NTU  =  —  (34) 

{-'min 

Hence,  the  heat  exchange  rate  between  the  hot  and  the 
cold  gas  streams  is: 


<Uv 


2F  . 

j .  v/Ohm  3”  h Act, a  4"  h Act,c  4~  *?C  one,  a 


4"  hConc.c) 


(28) 


(/  —  £(/max  (35) 

where  the  theoretical  maximum  heat  transfer  rate  is 


3.2.  Heat  exchanger  model 

3.2.1.  Design  of  heat  exchanger 

The  effectiveness-number  of  transfer  units  (e-NTU) 
method  is  used  to  model  the  pre-heaters  in  H2-  and  CH4- 
fed  fuel-cell  systems  [11].  In  this  study,  counter-flow  types 
of  heat  exchanger  are  selected  as  the  pre-heaters,  and  are 
sized  according  to  the  following  prescribed  conditions  given 
in  Table  2. 

3.2.2.  Performance  calculation 

The  simulation  code  is  designed  to  allow  the  users  to  vary 
the  inlet  and  operating  conditions  of  the  system.  To  deter¬ 
mine  the  outlet  temperatures  of  the  pre-heaters,  which  vary 
with  the  inlet  conditions,  the  heat  capacity  rates  of  the  cold 
and  hot  gas  streams  are  calculated: 

Cc  =  £  hccp_c  (29) 

Ch  =  £  "hCP4  (30) 

By  comparing  their  values,  the  lower  and  higher  values 
are  assigned  as  Cmin  and  Cmax,  respectively.  The  ratio  of  heat 


(/max  —  6'mjn  (7'h.i  Pq.\  ) 


(36) 


Based  on  the  energy  balance,  the  exit  temperatures  of  the 
hot  and  cold  gas  streams  from  the  heat  exchanger  are: 

Th.o  =  Th.i  -  X  (37) 

Tc.0  =  rc4  +  X  (38) 

'-'C 

Since  the  inlet  temperatures  of  the  heat  exchangers  are 
affected  by  the  exit  flow  temperature  of  the  SOFC  stack,  the 
above  procedures  are  repeated  until  the  exit  temperatures  of 
the  pre-heaters  converge. 


4.  Results  and  discussions 

Details  of  various  overpotentials  as  a  function  of  current 
density  in  a  fuel-cell  operated  at  800°C  are  shown  in  Fig.  4. 
The  limiting  current  density  at  this  operating  temperature  is 
2980  A  m  2.  Results  show  that  cathode  activation  and 
Ohmic  overpotentials  are  responsible  for  the  major  losses 
in  the  fuel-cell  over  the  normal  operating  range.  The  cathode 
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concentration  overpotential  becomes  significant  when  the 
fuel-cell  operates  near  to  the  limiting  current  density.  Com¬ 
pared  with  the  anode,  the  cathode  exhibits  higher  activation 
overpotential,  which  is  due  to  the  poor  ‘apparent'  exchange 
current  density  at  the  electrode/electrolyte  (LSM-YSZ/ 
YSZ)  interface,  as  will  be  seen  later  in  Fig.  8.  Since  the 
cathode  exchange  current  density  directly  affects  the 
electrochemical  reaction  rate  at  the  cathode,  it  can  be  under¬ 
stood  that  the  low  electrochemical  reaction  rate  in  the 
cathode  leads  to  high  cathode  activation  polarization  in 
the  fuel-cell. 

The  effect  of  temperature  on  the  resistance  of  the  anode, 
cathode,  electrolyte  and  interconnector  is  presented  in 
Fig.  5a-d,  respectively.  The  resistances  of  these  fuel-cell 
components  are  determined  by  the  resistivity  of  the  materi¬ 
als  used  and  their  respective  thickness.  The  results  show  that 
the  resistances  of  the  cathode,  electrolyte  and  interconnector 
decrease  with  increase  in  temperature.  By  contrast,  the 
anode  resistance  displays  the  opposite  trend.  The  resistances 
of  both  electrodes  are  negligible,  however,  compared  with 
those  of  the  electrolyte  and  the  interconnector,  despite  the 
fact  that  a  thick  electrode  (cathode)  was  used  as  the  fuel-cell 
support.  The  details  of  the  resistivity  and  thickness  of  each 
fuel-cell  component  are  listed  in  Table  3. 

The  performance  of  the  fuel-cell  as  a  function  of  tem¬ 
perature  is  presented  in  Fig.  6  in  the  form  of  cell  voltage 
versus  current  density.  Results  show  that  the  higher  the 
temperature,  the  lower  will  be  the  Nernst  potential.  As  more 
current  is  drawn  from  the  fuel-cell,  however,  a  higher  cell 
voltage  can  be  maintained  under  higher  temperature  opera¬ 
tion  at  the  same  current  density.  There  is  significant  reduc¬ 
tion  of  both  Ohmic  and  activation  losses  (due  to  increased 
exchange  current  densities)  at  high  temperature  operation. 


Table  3 


Resistivity  and  thickness  of  cell  components 


Material  used 

Ni-Y  SZ/Y  SZ/LSM-Y  SZ 

Anode  thickness  (|tm) 

150 

Anode  Ohmic  resistance  constant 

a  =  0.0000298%  =  -1392 

Cathode  thickness  (m) 

2E-3 

Cathode  Ohmic  resistance  constant 

a  =  0.0000811,  b  =  600 

Electrolyte  thickness  (|tm) 

40 

Electrolyte  Ohmic  resistance  constant 

a  =  0.0000294,  b  =  10350 

Interconnect  thickness  (|tm) 

100 

Interconnect  Ohmic  resistance  constant 

a  =  0.001256,  b  =  4690 

This  causes  an  improvement  in  fuel-cell  performance 
despite  the  lower  Nernst  potential  at  higher  temperature. 

The  effect  of  temperature  on  exchange  current  density  is 
shown  in  Fig.  7.  Both  the  anodic  and  cathodic  exchange 
current  densities  increase  exponentially  with  temperature, 
which  demonstrates  the  critical  dependency  of  this  perfor¬ 
mance  parameter  on  temperature.  The  results  also  clearly 
indicate  that  the  exchange  current  density  of  the  anode  is 
more  sensitive  to  the  operating  temperature  than  that  of  the 
cathode;  the  difference  between  the  two  current  exchange 
densities  is  about  62%  at  1073  K. 

The  performance  of  the  fuel-cell  as  a  function  of  operating 
temperature  is  presented  in  Fig.  8  in  the  form  of  power 
density  versus  current  density  traces.  Results  show  that  the 
higher  the  temperature,  the  higher  is  the  peak  power.  The 
location  of  the  peak  power  shifts  towards  higher  current 
density  as  the  temperature  increases,  which  shows  an 
improvement  in  limiting  current  density  at  higher  tempera¬ 
ture.  As  a  compromise  between  running  a  fuel-cell  at  peak 
power  and  the  stability  of  operation,  it  is  advisable  to  operate 
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Fig.  5.  Effect  of  temperature  on  (a)  anode  Ohmic  resistance;  (b)  cathode  Ohmic  resistance;  (c)  electrolyte  Ohmic  resistance;  and  (d)  interconnector  Ohmic 
resistance. 


the  SOFC  at  a  current  density  slightly  less  than  that  corre¬ 
sponding  to  the  peak  power. 

The  performance  of  the  fuel-cell,  operated  between  1073 
and  1273  K,  is  presented  in  Fig.  9  in  the  form  of  maximum 


power  density  versus  the  percentage  of  current  density 
corresponding  to  the  peak  power  at  the  respective  limiting 
current  density.  Note  that  the  ratio  of  current  density  to 
limiting  current  density  corresponds  to  the  fuel  utilization 
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Fig.  6.  Effect  of  temperature  of  the  cell  voltage  vs.  current  density  trace. 
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Fig.  8.  Effect  of  temperature  on  power  density. 


Fig.  9.  Peak  power  density  output  at  respective  fuel  utilization  for  temperatures  between  1023  and  1273  K. 
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Fig.  10.  Effect  of  molar  water/fuel  ratio  on  heat  required  for  reforming  (fuel  utilization  rate:  75%,  reformer  efficiency:  90%,  973  K  (fuel),  and  873  K  (air)). 


rate.  The  trend  for  the  maximum  power  is  expected  to 
increase  with  increase  in  operating  temperature,  but  such 
a  trend  is  limited  by  the  material  constraints  and  the  con¬ 
centration  polarization.  As  such,  it  is  impossible  for  the  fuel¬ 
cell  to  run  at  a  fuel  utilization  rate  of  100%. 

The  effect  of  the  amount  of  water  used  in  the  steam 
reforming  on  the  heat  requirement  is  shown  in  Fig.  10.  The 
general  trend  is  that  the  heat  required  to  reform  CH4  gas  into 
H2  through  steam  reforming  increases  with  increase  in  the 
water/fuel  ratio.  As  mentioned  earlier,  a  molar  water/fuel 
ratio  of  2  is  sufficient  to  reform  the  CFt4,  but  more  water  is 
normally  used  in  order  to  prevent  undesirable  formation  of 
carbon  due  to  possible  thermal  decomposition  of  the 
methane.  In  addition,  since  the  current-reforming  model 
is  based  only  on  fixed  chemical  composition  of  the  product 
gases  and  is  independent  of  temperature,  no  further  insight 
can  be  gained  from  this  result. 

The  enthalpy  and  exergy  at  each  state  point  of  the  H2-fed 
SOFC  system  are  presented  in  Fig.  11.  In  this  particular 
study,  the  SOFC  stack  operates  at  a  75%  fuel  utilization  rate 
and  at  1  atm  for  both  the  anode  and  cathode  flow.  The  fuel 
and  air  inlet  temperatures  of  the  stack  are  set  to  973  and 
873  K,  respectively.  The  performance-related  data  enthalpy 
and  exergy  are,  respectively,  normalized  by  the  LHV  and 
chemical  exergy  of  H2.  In  addition  to  these  data,  the 
temperature  at  each  state  point  is  also  shown  in  Fig.  1 1 . 
The  first  value  in  each  control  volume  is  the  exergy  destruc¬ 
tion  while  the  second  value  is  the  exergy  destruction  normal¬ 
ized  by  the  chemical  exergy  of  H2.  Energy  transfer  through 
the  system  boundary  of  the  SOFC  stack  includes  both  heat 
and  work  transfer.  The  former,  as  mentioned  earlier,  is  due 
mainly  to  the  irreversible  electrochemical  reaction,  while  the 
latter  is  the  output  electrical  energy  of  the  SOFC  stack.  Both 
heat  transfer  and  work  transfer  are  also  normalized  by  the 
LHV  and  chemical  exergy  of  H2,  respectively.  If  a  system 
boundary  is  drawn  to  enclose  all  the  control  volumes  and 
allow  the  waste  flow  energy  and  heat  transfer  discharge/ 


dissipate  to  the  environment,  it  is  obvious  that  the  summa¬ 
tion  of  inlet  and  outlet  enthalpies  should  obey  the  conserva¬ 
tion  of  energy.  In  the  case  of  exergy,  the  input  exergies 
should  be  balanced  by  the  outlet  exergies  plus  the  exergy 
destructed  in  all  control  volumes.  Hence,  summation  of 
enthalpies  at  state  points  1  and  2  should  be  equal  to 
enthalpies  at  state  8  plus  the  heat  dissipated  to  the  environ¬ 
ment  by  the  afterburner  and  SOFC  stack  and  the  electrical 
energy  output  from  the  system.  Likewise,  the  summation  of 
exergies  at  state  points  1  and  2  should  be  equal  to  the  exergy 
at  state  8  plus  exergy  loss  due  to  heat  transfer  from  the 
afterburner  and  SOFC  stack,  electrical  exergy  output  from 
the  system  and  the  exergies  destructed  in  all  control 
volumes.  In  the  SOFC  stack,  the  normalized  values  of  the 
electrical  work  are,  respectively,  the  first  law  and  second 
efficiencies.  In  the  case  of  a  H2-fed  system,  the  second  law 
efficiency  is  higher  than  the  first  law  efficiency  due  to  the 
fact  that  the  LHV  of  H2  is  higher  than  its  corresponding 
chemical  exergy.  A  study  of  the  system  shows  that  the 
highest  energy  lost  is  due  to  the  heat  transfer  at  the  after¬ 
burner,  which  contributes  25.75  and  19.06%  of  the  LHV  and 
chemical  exergy  of  H2,  respectively.  It  can  be  explained  that, 
for  system  of  this  kind,  the  isothermal  afterburner  (i.e. 
maximum  heat  loss)  is  used  to  prevent  excessive  increase 
in  the  temperatures  of  H2  and  air  at  pre-heaters  2  and  1, 
respectively,  so  that  the  SOFC  stack  can  operate  within  the 
allowable  feedstock  temperatures.  Note,  with  too  high  a 
feedstock  temperature,  there  will  be  an  adverse  effect  on  the 
reliability  of  the  stack  operation  if  the  stack  temperature 
operates  far  from  its  design  condition.  To  improve  the 
overall  system  efficiency,  co-generation  may  be  applied  to 
make  use  of  the  waste  heat,  which  is  otherwise  lost  through 
the  cooling  system  to  maintain  the  isothermal  condition  of 
the  afterburner.  In  view  of  the  low  exergy  destruction  in  all 
control  volumes,  an  improvement  of  thermal  efficiency  can 
be  obtained  by  operating  the  stack  at  high  fuel  utilization, 
but  not  excessively  high  to  cause  problems  of  severe 
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Fig.  11.  Hydrogen-fed  system  operating  at  75%  fuel  utilization,  1  atm  (anode  and  cathode),  973  K  (fuel),  and  873  K  (air). 


concentration  overpotential  and  overheating  at  high  stack 
temperatures. 

Slightly  more  complicated  system  is  shown  in  Fig.  12  and 
consists  of  a  vaporizer,  a  mixer  and  a  reformer  in  addition  to 
the  basic  system  requirements  described  above.  To  provide  a 


comparison  with  a  FF-fed  system,  all  performance-related 
input  parameters  are  set  to  the  same  values.  Results  show 
that  the  first  and  second  law  efficiencies  of  this  CH4-fed 
SOFC  system  are  62.19  and  59.96%,  respectively,  which  are 
higher  than  those  in  a  H2-fed  system.  Note,  the  first  law 


T1 1:552. 37K 


Fig.  12.  Methane-fed  system  operating  at  75%  fuel  utilization,  1  atm  (anode  and  cathode),  973  K  (fuel),  and  873  K  (air). 
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efficiency,  in  this  case,  is  higher  than  the  second  law 
efficiency  and  is  attributed  to  the  fact  that  the  chemical 
exergy  of  the  fuel  is  higher  than  its  corresponding  LHV.  In 
this  system,  all  losses  (exergy  destruction  in  all  control 
volumes,  heat  transfer  and  flow  exergy  leaving  the  system) 
are  rather  well  distributed  with  higher  flow  exergy  loss  at  the 
exit  of  pre-heater  1  (9.57%)  and  exergy  destruction  in  the 
vaporizer  (7.22%).  Compared  with  a  H2-fed  system,  the  heat 
transfer  loss  to  the  environment  in  a  CH4-fed  system  is 
lower.  This  is  attributed  to  a  higher  demand  for  heat  energy 
used  by  the  reformer  and  vaporizer. 

5.  Conclusions 

A  simulation  program  has  been  developed  to  simulate  two 
simple  SOFC  power  systems  which  use  hydrogen  and 
methane  as  the  feedstock  fuels.  Based  on  this  study,  the 
following  conclusions  are  made. 

•  The  fuel-cell  model  developed  in  this  study  can  replicate 
well  the  characteristics  of  the  fuel-cell  under  different 
operating  conditions,  as  well  as  and  the  electrical  proper¬ 
ties  of  the  materials  used  in  the  cell  components.  The 
parameters  related  to  electrochemical  processes  and 
transport  phenomena  are  all  expressed  as  the  function 
of  thermodynamic  state. 

•  The  fuel-cell  model  can  be  used  to  simulate  different 
designs  of  SOFC  provided  that  the  geometric  specifica¬ 
tions,  performance-related  parameters  and  operating  con¬ 
ditions  are  specified. 

•  In  a  simple  SOFC  power  system  with  only  waste  heat 
recovery  used  for  pre-heating  the  fuel  and  air,  the  system 
efficiency  can  be  improved  by  operating  the  fuel-cell 
stack  at  a  high  fuel  utilization  rate,  but  not  excessively 


high  to  cause  problems  of  concentration  overpotential  or 
overheated  cells  due  to  high  stack  temperature  associated 
with  high  cell  polarization. 

•  At  a  75%  fuel  utilization  rate,  1  atm  of  anode  and  cathode 
flow  pressure,  a  973  K  of  fuel  inlet  temperature  of  973  K 
and  an  air  inlet  temperature  of  873  K,  simulated  results 
showed  that  a  FF-fed  system  can  achieve  50.97  and 
52.28%  of  the  first  and  second  law  efficiencies,  respec¬ 
tively.  For  a  CH4-fed  system  under  the  same  settings,  the 
corresponding  first  and  second  law  efficiencies  are, 
respectively,  62.19  and  59.96%. 
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